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AXIALLY SYMMETRIC SHAPES WITH MINIMUM WAVE DRAG 1 

By Max. A. HEASLETand Feanklyx B. Fcllee 


SUMMARY 

The external wave drag of bodies of revolution moving at 
supersonic speeds can be expressed either in terms of the geometry 
of the body, or in terms of the body-simulating axial source 
distribution. For purposes of deriving optimum bodies under 
various given conditions, it is found that the second of the 
methods mentioned is the more tractable. By use of a quasi- 
cylindrical theory, that is, the boundary conditions are applied 
on the surface of a cylinder rather than on the body itself, the 
variational problems of the optimum bodies having prescribed 
volume or caliber are solved. The streamwise variations of 
cross-sectional area and drags of the bodies are exhibited, and 
some numerical results are given. The solutions are found to 
depend upon a single parameter involving Mach number and 
the radius-length ratio of the given cylinder. Variation of this 
parameter from zero to infinity gives the spectrum of optimum 
bodies, for the prescribed condition, from the slender-body result 
to the two-dimensional. The numerical results show that for 
increasing values of the parameter, the optimum shapes quickly 
approach the two-dimensional. 

A reciprocity relation for axial flow is derived, and it is used 
in formulating the variational problems in terms of the drag 
formula involving geometry. Formulation of the minimum 
problems in terms of combined flow fields is found to lead to 
extremely simple relations that are satisfied by the flow field 
induced by optimum bodies. The combined flow concepts are 
also useful, for example, in checking results found by other 
means. 

INTRODUCTION 

The design of minimum-drag configurations is one of the 
fundamental problems of aerodynamics. For many en- 
gineering purposes it is, furthermore, possible to make useful 
predictions and design calculations for steady flight by con- 
sidering additively the drag attributable to the viscous 
nature of the air and the drag that occurs in an inviscid 
medium. Since efficient flight is closely associated with the 
use of aerodynamic shapes producing relatively small dis- 
turbances in the air, the analysis upon which the inviscid- 
fluid theory is based can, in many cases of practical interest, 
be further limited to first-order approximations involving 
small perturbations. For supersonic flight speeds such an 
analysis is linear, the perturbation velocity potential of the 
flow field satisfies the wave equation, and the pressure drag 
of nonlifting configurations results from the accumulation of 
energy in the waves induced by the body during its motion. 

The purpose of the present paper is to show how most 


favorable body shapes, under various given conditions, can 
be derived by using formulae for drag prediction that are 
based upon the linearized theory. The type of body to be 
treated is a nacelle- or duct-like configuration (nonlifting and 
having axial symmetry) which induces perturbations that 
are specified on the surface of a circular cylinder. The 
analysis might be termed quasi-cylindrical, since boundary 
conditions are applied on the surface of a cylinder rather 
than on the body itself. Only the external flow is considered. 

There are two rather different methods available for the 
calculation of drag of such bodies. The first, given by Ward 
in reference 1, expresses the drag in terms of the geometry of 
the body and of a weighting function first encountered by 
Lighthiil (ref. 2) in connection with the drag of fusiform 
bodies. The second result, published recently by Parker 
(ref. 3), is a formula in which the drag. is expressed in terms 
of the strength of an axial source distribution that simulates 
the body shape. Generally speaking, the formula giving 
drag directly in terms of geometrical characteristics would be 
preferable, since the usual auxiliary conditions in variation 
problems, such as given volume, given caliber, etc., are also 
expressed in geometrical terms. Unfortunately, however, 
the variational problem in this case leads to an integral 
equation whose kernel is the Lighthiil function mentioned 
previously, and the properties of this function are not at 
present well enough known to enable one to solve the integral 
equation by other than numerical methods. On the other 
hand, the expression for drag in terms of sources leads to a 
tractable integral equation, although the relations between 
source strength and geometry are somewhat complex. 

Problems of the sort to be treated here have been attacked 
by Ferrari (refs. 4 and 5) and by Parker (ref. 3). The first- 
named author has approached the problem of minimum drag 
with assorted isoperimetric conditions by both the above- 
mentioned methods, but the main effort was made in con- 
nection with the source-strength method applied in con- 
junction with a control surface consisting of a frustum of a 
cone. A large number of cases have been worked out, 
mostly by numerical methods. The other work, reference 3, 
gives a solution to the problem of the 'minimum-drag body 
with given caliber, making use of boundary conditions on the 
Stokes’ stream function, rather than the potential function. 

In this paper we shall approach the problem by the use of 
both methods outlined above. In an introductory section, 
the operational approach to the wave equation is extended 
to bodies having peripheral as well as longitudinal variations 


1 Supercedes NA O A TN 3389 by Mai. A. Hcaslet and Franklyn B. Fuller, 1955. 
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of surface shape. The analysis is then restricted to the case 
of axial symmetry and the two drag formulae are given. 
Then a reciprocity relation for axial flow is derived, and the 
notion of combined flow fields is introduced. This device 
leads, through application of the reciprocity relation and the 
drag formula in terms of body geometry, to extremely simple 
physical characterizations of the flow fields associated with 
optimum bodies. Next, in order to derive explicit expres- 
sions for some optimum bodies we consider the source- 
function approach in combination with a cylindrical control 
surface on which boundary conditions are specified. The 
results obtained are discussed with the aid of numerical 
examples, and, finally, the reciprocity relations derived 
earlier are exhibited in terms of the explicit solutions found, 
and some uses of the reciprocity results are indicated. 

The appendix is devoted to summarizing the results of the 
minimizations for the convenience of the reader. 
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LIST OF IMPORTANT SYMBOLS 


speed of sound in free stream 
strength of source distribution 
function used in isoperimetric problems (See 
eqs. (60).) 


pressure coefficient, 

2 ° 

drag 

complete elliptic integral of second kind of 
modulus k 

modulus of elliptic integrals 
complete elliptic integral offirst kind ofmodulusi 
Bessel functions of order m (See ref. 8.) 
length of body ~ 


Mach number in the free stream, 


Uo 

(Lq 


direction cosines with respect to Cartesian axes 
of the inward normal to a surface 
pressure 

pressure in the free stream 
pressure in a combined flow field, y—y 

dynamic pressure, „ pJJ 3 

radial coordinate, VP+z 1 

incremental radius on control cylinder due to 
source distribution along axis 
radius of cylindrical control surface 
cross-sectional area of a body 
S(x)-S(o) 

Cartesian coordinates 

perturbation velocities in x,y,z directions, respec- 
tively 

free-stream velocity 

perturbation velocity in radial direction 
volume of body 

additional volume wrapped on cylindrical con- 
trol surface. . 

function defined in equation (25) 
parameter of elliptic integral of third kind 
Mo 3 — 1 
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dimensionless streamwise coordinate, j 
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angular coordinate, tan -1 - 

Lagrange multipliers in isoperimetric problems 
complete elliptic integral of third kind of modu- 
lus k and parameter a 3 (in notation of ref. 20) 
free-stream density 
§R 
l 

perturbation velocity potential 

SUFFIXES 


differentiation with respect to streamwise coordi- 
nate 

quantity evaluated in reversed flow field 

ole. 


Laplace transform 

V 

dimensionless quantity as V*=-^> 


INTRODUCTORY ANALYSIS 

The analysis to be given here is adapted to boundary con- 
ditions specified on a right circular cylinder so oriented that 
its axis is parallel to the free-stream velocity vector. 
Immediate application thus follows for quasi-cylindrical 
shapes that deviate slightly, both longitudinally and periph- 
erally, from a cylindrical control surface although the expres- 
sion for drag can be extended to include the domain of 
slender-body theory. 

Consider a fixed Cartesian coordinate system in a super- 
sonic free-stream of velocity U 0 and Mach number M 0 — U„/a 0 
>1 where a 0 is the velocity of sound in the free stream . The 
x axis is alined with the direction of the flow and the lateral 
coordinates y,z may a lso be expressed in polar coordinates 
r,6 where r=Vy*+ 2 *, 0—tan~ l z/y. A cylindrical control 
surface of radius r=/i=const. is given with the range 
0 <x <1 and on this control surface the perturbation velocity 
components, together with their gradients, are small relative 
to U 0 and U 0 /l. Under these conditions the field external 
to the cylinder of radius R has for its governing equation the 
linear relation 

/SV«— <p n —Va= 0 (1) 

where the subscript notation denotes partial differentiation, 
<p{x,y,z) is the perturbation velocity potential yielding the 
perturbation velocity components 

u(x,y,z)=<p x (x,y,z), v(x,y,z)=<p s (x,y,z), w(x,y,z)=<p t (x,y,z) 

and jS 3 =M„ 3 — 1. The boundary conditions on the body aro 
to be taken in the form 


<Pr{x,r,e)\ TaR =U 0 6(xfi), 0<x<l (2) 


where 0 is a known function of x and 6. 


A GENERAL SOLUTION OF THE WAVE EQUATION IN CYLINDRICAL COOR- 
DINATES 

If equation (1) is rewritten in the form 


0*Vxx— <Prr— (l/r)Vof>=0 (3) 

it is possible, through separation of variables, to derive a 
general solution representing a rectilinear distribution o 
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source and multipole singularities. Tliis general solution can 
be found by use of tbe Laplace transformation. By defini- 
tion, the Laplace transform 2 (see ref. 6) of a function 
F(x,r,d) is F(s',r,6) where 

n cd 

F(s]rfi)= \ e~‘ I F(x,r,6)dx (4) 

If one employs this transformation and applies initial condi- 
tions consistent with supersonic flow theory (ref. 7), equation 
(3) becomes 


^^~<p rr ~(l/r)<p r —(lfry<pee=0 (5) 

The tranform of the perturbation velocity, potential is 
assumed separable in the form 

p(s;r,0)=f(r,8) cos md 

and it follows directly that t(r,s) must satisfy the ordinary 
differential equation 

_k ri+— i r=o 

d(0rs)^(3rs d(J3rs) [_ G9™) 5 J * U 

Thus, the solution can be written 


£(s;r,0)=— A-; jp COS md [A m (s) K n (prs)+B m (s) I n (firs)] 


where K„ and I m are modified Bessel functions in the nota- 
tion of reference 8. The asymptotic expansions for the 
Bessel functions show that I m yields incoming waves suit- 
able for the analysis of flow inside a tube or cylindrical 
control surface; K m yields outgoing waves that are suited to 
the calculation of the field external to a tube. It follows 
that one has, in the latter case, 

v(s;r,0)=— ^-J^A„(s)E m (fir8)cosme (6) 

0 

The inversion of equation (6) can be achieved in two 
ways. First, from reference 9, page 277, and the convolu- 
tion integral, one gets 


¥>(r,r,0)=- 


2 TT 
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~ pr A n (xi)dxi 


o ^(x-xtf—pr 2 


r *-pr A <n( x i) cosh ( m cosh" 1 ) dXi 

s cos md \ , v ^ 

i Jo Jfx—XiY—B^r 2 ■ 


V(Z-X l) 2 -^ 
Second (see, e. g., ref. 8, p. 79), one has 

K,m 

Thus equation (6) can be rewritten as 
W,r,6)=-^ jb^(s).£: o (fire)+ 


£(-j)”cos nU>(£)‘^KOr.; 


(7) 


and the inversion is 


<p(x. 




~ <3r Ajx^Axi 


where the function C m (x) is given by (from operational 
calculus rules) 

C m (x)= f cfo, • r • \ dx* \ A„ fo) d x i 

Jo Jo Jo 


(m 


J o (x—x^-'AJx^dxi 


0) 


Equation (8) expresses the solution in the usual form, 
given, for example, in reference 10, page 527. For some 
purposes, numerical calculations for example, equation (7) 
has advantages over equation (8). The two solutions 
express the perturbation velocity potential in terms of 
distributions of singularities along the central axis, the first 
term representing a distribution of supersonic sources of 
strength A 0 (x)dx, and the subsequent terms representing 
multipoles of order m. 

It is of interest to calculate the limiting forms of equations 
(7) and (8) for large and small values of r. For large r, 
equation (6) becomes 


vis;r, 0) ~ jp A m {s)^- e~ Br ‘ cos md (10) 

where the asymptotic form 


K m {z) 



e 1 


has been used. The perturbation velocity potential is then 


<p(x,r,d)' 


— y) cos md 

2-ry2/3r o 


Jo " 


ri„(S) 


dS 


-Jx—t—pr 


(ID 


The ultimate attenuation of <p with lateral distance is there- 
fore fixed by the factor 1 /-y/r. For small r, equation (6) be- 
comes 


7p(8-,r,d)~—^—A 0 {8) (in ^+t)+ 

2^!)( m -i)i(A)"cos 7 * 0 ] (12) 

where 7=0.577 is Euler’s constant. The inversion of equa- 
tion (12) is 

<p{x,r,d)~ ^\A 0 (x)ln I A 0 ^ln\x—Xi\dxi+ 

±, (|)" 003 ^ ^ (X) ] (13) 

where O m { x) is defined in equation (9). This result was used 
by Ward (ref. 11) as a basis for the development of slender- 
body theory. 


* It will be assumed through the present section that the origin lies upstream of all disturbance points In tbe flow field. Subsequently, the origin will be shifted so as to lie at the upstream 
face of the control surface or body. 
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As presented, the above general solutions (eqs. (7) and (8)) 
were not related to specific boundary conditions. The formal 
development of this relation is straightforward and leads to 
an explicit solution for boundaiy conditions given on the 
cylindrical control surface at r=l?=const. Let the given 
conditions be 


<Pt]t-r=U 0 G(x,B)=U b S gm(x) cos md 
0 


(14) 


From equations (6) and (14), one has 
Vr(8',r,0)]r-B=— ^ ib A„(«) X^OSrs)^ ^ cos me 

co 

=U 0 2D In.(«) cos me (15) 

0 

Since 

^ K„(firs)=fisK wl ' (firs) 
equation (15) yields 

iw-f '; wm <16) 

and the transformed velocity potential is, from equation (6), 


U 0 ^ g m (s) K m (J}rs) 




(fiRs) 


cos md 


(17) 


In order to give the desired expression for <p(x,r,6 ) it is neces- 
sary to calculate the inverse Laplace transform of the func- 
tions K m (firs)/K m '(fiRs). This task has been undertaken by 
Mersman (ref. 12). 

EXTERNAL WAVE DRAG OP QUASI-CYLINDRICAL BODY OP REVOLUTION 
IN TERMS OF ITS GEOMETRY OR SOURCE DISTRIBUTION 

Attention is now restricted to flow fields possessing axial 
symmetry with respect to the steam direction. Independence 
with respect to 0 then reduces equations (7) and (8) to 


** ’ 27tJo V(x— x,) 2 — /3¥* 


(18) 


and the velocity potential is expressed as a rectilinear dis- 
tribution of supersonic source potentials. Operationally, 
equation (18) takes the form 


<p(s;r)= ~~ A p (8)K 0 ($r8) 


(19) 


The axes may now be considered as shifted so that the source 
distribution starts at x=—f3R and induces perturbation 
velocities on the cylindrical surface r—R, OAx^l. For 
r>R one then has the disturbance field associated with a 
body of revolution that deviates only slightly from the 
cylinder r—R. The wave drag of such a body can then be 
expressed in two ways: first, as a function of the body geom- 
etry; second, as a function of the source-strength distribu- 
tion. The first result has been given in reference 1. To the 
order of accuracy to which this control-surface theory applies, 
the slope of the resulting surface is 


dr 

dx' 


,STx) 
2ttR ' 


1 " 
■-fjVr 

U o r- 


( 20 ) 


where S'(x) is the stream wise derivative of local cross- 
sectional area of the body. This condition, together with 
equation (19), yields 


A 0 (s)_ S'(s) 

U 0 QRsK^ifiRs) 


where S' («) means the Laplace transform of S' (x), and 
yjs-R) S'(8) K o (0Rs) 


U 0 


2irftRs KiifiRs) 


( 21 ) 


( 22 ) 


In order to calculate drag, pressure on the body is next 
evaluated. Denoting by p and p„ local and free-stream 

pressure and setting $r 0 =| poU 0 2 , one has in linearized theory 


<Jo Jr- 1 


2 u(x,R) 


U 0 


(23) 


From equation (22) 
u(s;R) 


U„ 


1 cr/r x KpifiRs) 

2irfiR ° w KMis) 


= _ S'{s ) r Ki(pRs)-K 0 (pR8) J , . 

2*m L KtfRs) J { v 

The inveise transform of the second term involving the 
Bessel function leads to the function W{x) introduced by 
Ligh thill (ref. 2). By definition, its transform is 


W(s) 

K&) 


(25) 


Pressure distribution on the body can then be calculated 
from the expression 

V- [ s ' ( *>-r ' s ' w Kiir) ml m 

The function W(x) is shown in figure 1; tabular values for 
— 2<x<10 are given in reference 1. 
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The external wave drag C Do of the body is finally deter- 
mined by direct integration 


n drag 1_ C l P—P„ 

D ° irR 1 ^ o tR 2 J 0 q 0 


S'(x)dx 


(27) 


and from equation (26) is 


In a later section entitled “Geometric Criteria for Minimum 
Drag,” the role equation (28) plays in problems involving 
drag minimization will be discussed. For the present, it 
may be remarked that although the magnitude of the in- 
fluence function W(x) is known, its analytic properties are 
not well enough defined to permit easy manipulation. It 
will become more apparent later that for certain minimum- 
drag problems an advantage is provided when one deals 
directly with source distributions and establishes the rela- 
tionship between geometry and source strengths as a 
separate part of the analysis. 

Equation (18) expresses the potential of a source distribu- 
tion of strength A 0 (x). On the cylindrical control surface 
r=R and within the range 0 <x<l an effective body shape 
is induced and the drag of this body can be calculated as 
follows. The stream wiso and lateral perturbation-velocity 
components are, respectively, 


, v 1 

C x -p A/ix^dXi 

(29) 

«*<*’’•> — & 

J-/s« aJix—Xx y—pr* 


(x—x 1 )A < ,'(x,)dx l 
-pa ■\,/(x—X 1 ) s — jSV 1 

(30) 


where A o (x l )=0 for Xi <—f}R. The effective body, within 
the range 0 <x<l, is fixed by the boundary conditions of 
equation (20) and its external wave drag is 

D =— 2ttp,,rJ <p z (x, R)<p,(x,R)dx 

Pa C 7 r x -^ A o'( x 0( x — x i)<foi C x ~^ AJjx^dx^ fn ^ 

~2vJo v(*-*i y-FR 2 J -PB m 


The dummy variables X\, x* can be interchanged; if one then 
combines the two expressions of equation (31) and inverts 
the order of integration, 3 the integration with respect to x 
can be performed and there results 


x 


PR 
l- PR 


Ao\x^dx i 


AJ(x, j) cosir 


(i-xoa-xa-pip 


f}R(xi—Xi) 


dXi (32) 


as given in reference 3. 

It is of interest to remark that although equation (32) 
uses only a knowledge of the function A 0 (x) in the range 
—fiR<x<l—@R, the drag that is calculated presupposes a 


specific source distribution function in the range l — fiR<^x 
if one wishes to identify the drag with a geometric shape. 
Thus, as in figure 2, if the body shape near r=R is assumed 
to have some arbitrary variation for 0 <x<l, and to 
straighten out into a purely cylindrical surface downstream 
of x=l, a source distribution function is required downstream 
of x=l—pR to produce the cylinder. 


✓ 

/ 



Figure 2. — Body induced on a control surface by an axial source 

distribution. 


The fact that the stream velocity is supersonic means 
that upstream influences of A 0 (x) for x>l—@R cannot be 
felt on the body and explains why the drag of a complete 
geometric shape can be determined from its source distribu- 
tion without knowing the complete details of the distribution 
function. 

As another example of the use of equation (32) consider, 
as in figure 3, a circular body extending from x=—$R to 
x=l with a cylindrical afterbody of radius R aft of x=l. 
If the source distribution of this body is known as, say, for 
example, in the case of a cone or slender body of revolution, 
the body drag can be determined by using tho surface 
r=R,0<x<l as a control surface and calculating momen- 


r 



Figure 3. — Body nose induced by an axial source distribution. 


* Tho Inversion of order is permissible only if A»(x) is suitably well behaved, a point that will arise later in the determination of the optimum body with given caliber. 
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turn transport through the control surface. Equation (32) 
is the exact expression for the body drag, and, again, requires 
no knowledge of source strength beyond x=l—fiR. 

COMBINED FLOW FIELDS 

One method of attack that has proved to be extremely 
helpful in the analysis of problems in aerodynamic theory 
involves a symmetrization process in which flow fields in 
both forward and reverse flow are related. Attention, up 
to the present time, has been devoted principally to planar- 
type problems and in reference 13 Jones has used this 
approach to derive criteria that appear in the minimization 
of wave drag of, for example, nonlifting wings having speci- 
fied thickness ratios or volumes. In this section, a brief 
discussion is given, using the methods of reference 14, of the 
way these concepts appear in cylindrical-control-surface 
analysis. 


If* li+"> sf ^ If) * 

= IIK - ‘^s+^lj; + ^§;)' !z (34) 

applies where the surface integration extends over Sj, 2 2 , 2 3 
and ni, 7i2, n s are direction cosines, with respect to the x, y, z 
axes, of the surface normal directed inward into the region. 

It is customary to re-express relations like equation (34) 
in terms of a newly defined directional derivative along a lino 
termed the conormal. In this manner, the equation becomes 

JKsfHJO (36) 

where 

< 3 « 


THE RECIPROCITY RELATION FOR AXIAL FLOW 

Equation (1) can be written 


IAfp) s ^<Pxx~<Pn~ <Pi *— 0 (33) 

where L(<p) is a self-adjoint linear operator. Let now 
ip(x,y,z) and Zl(x,y,z) be two solutions of equation (33) 
satisfying boundary conditions given on a circular cylinder. 
Reciprocal relations betwen ^ and Q can be derived by 
applying Green’s theorem over a prescribed geometric 
region. Consider, as shown in figure 4, the cylindrical con- 
trol surface extending from x=0 to x—l and draw the 
enveloping Mach cones at the front and rear of the surface. 



Figure 4. — Surfaces of integration for combined-flow analysis. 


and the direction cosines v u r 2 , v z of the conormal are derived 
from 

— 71 1 /3 2 =Afi, ti 2 =A»' 2 , n z ~Av z 

By calculation of the respective normals n lt n 2 , n 3 and using 
the relation vf-\- j^ 2 +f 3 2 = 1 , it is readily found from the 
equations defining the conormal that on the surface 2j, the 
conormal is normal to the surface andA=l; on a Mach cone, 
the conormal lies along the cone and A=/S 

Let now ^ be set equal to <p(x, r, 0), the perturbation 
velocity potential associated with boundary conditions in a 
forward-flowing stream, and lot fi be u(x, r, 0), the x-wisc 
component of perturbation velocity associated with boundary 
conditions in a stream flowing in the reverse direction. 
Under' these conditions, equation (35) becomes 

I.” *•/.' * If ix+fl !! u If If* If 

-jr m £ * ? ■ & « if* i -'Wi- % is ‘ 

On the Mach cone 22 , the perturbation potential may arbi- 
trarily be set equal to zero and its conormal derivative along 
the cone will also be zero; as a consequence, the second terms 
on both sides of the equation vanish. Since the flow fields 
are irrotational, &u/dr=2)p r /dx where v, is radial velocity. 
After making this substitution and integrating the first term 
in the right member by parts, one gets 


R f dd f uVrdx=R f \ <p(l,R,6)v T {l,R,0)— f v T udx\dO~ 
Jo Jo Jo L Jo J 


The last integral becomes 


Denote the cylindrical surface by 2i, the front Mach cone 
x— (3 r=— by 2 2 , and the rear cone x+/S r=l+0R by 
2 3 . These surfaces enclose a toroidal region, bounded 
internally by 2i and externally by 2 a and 2 S . It follows 
from Green’s theorem that the integral relation 


«■*(!)* 


and for the given boundary conditions it is possible to show 
that along a conormal of 2 3 the relation t> r =/3 u holds and 
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•/r u is independent of v. The integral then can be rewritten 
os 


uf (£) dv=(3R£ r u(l,R,d)'p(l,R > e)d6 

and one has, finally, the desired reciprocal theorem 

—R f dO f u(x,R,d)v,(x,R,d)dx=R f do f u(x,R,d)v r (x,R,6)dx 
Jo Jo Jo Jo 

(37) 

It is not the purpose here to exploit the various applica- 
tions of equation (37) ; rather, the role played by the recip- 
rocal relation in drag calculations will be considered. In 
the forward and reverse flow fields, the pressure-velocity 
relations of linearized theory are 

■p—Po=—pJU 0 u, p—p 0 =p 0 U 0 u (38) 


The proof follows from the relations 

S'(x)=S'(x) 
u(x,R') 


Uo 

2tt(3R 

u(x,K) 

1 

U 0 

2w$R 

P(x,R) 

; 1 r,, 




^ ( 42 ) 


GEOMETRIC CRITERIA FOR MINIMUM DRAG 

Consider now the problem of minimizing the wave drag of 
a quasi-cylindrical body subject to the condition that the 
volume of the body is constant. The body surface may be 
defined by 


If, furthermore, thickness distributions of the form 
r=f(x,6), r=f(x,6) 

are placed on the cylinder r=R, the boundary conditions are 

1 agl _df _1_ d$fl <5/ 

U„ drJ r „ R — bx U 0 dr} T-R — dr 


r=f(x,d)=R+g(x,6) (43) 

The function g(x,6) determines the magnitude of the surface 
displacement from the cylinder r—R; these displacements, as 
well as their gradients, are assumed small and we also assume 

g(x,0)^O for riSO and l^x 


Equation (37) can then bo written 

-r **£<*-*> ■*/>-»> i* 

An immediate consequence of this last result is that for a 
unique thickness distribution, that is, for /(r,0)=/(r,0), 
the drag of a body in forward and reverse flow is the same. 
This follows from the fact that for quasi-cylindrical bodies 
the relations for drag are, respectively, 

For fixed geometry, therefore, drag is equal to half the sum 
of these two expressions 

Defining pressure P(x,r,6) in the combined flow fields by the 
following 


P{x,r,6)=p- 

-p=-pJJo(u+U) 

(40) 



(41) 


If the body has axial symmetry, equation (41) reduces to 
the form given in equation (28). To show this, one notes 
first that P and/ are independent of 6 and that equation (41) 
becomes 

D=\ ^PS'(x)dx 


If equation (41) is integrated by parts with respect to x, the 
wave drag of the body becomes 

D= : ^j\l0j l o P'(x,R ; O)g(x > e)dx (44) 


where the prime indicates x-wise differentiation of P. The 
imposed geometric constraint on the variational problem is 


\ f dB f f \x,B)dx^\ f dO f [R 2 +2Rg(x,d)] dx 
6 J 0 Jo 4J 0 Jo 

=wRH+R f dB f g(x,6) dx=F== const. (45) 
Jo Jo 


where V is the total volume. The problem thus becomes 
one of minimizing the expression 


D-pV=-R £dej l o P'g(x,d) dxp 

n JjrjRZ-f dQ^g(x,0) (46) 

where n is the Lagrangian multiplier. Carrying out the 
variation, one has 


S(D-!iV)=~£ T dsJjP'^+giSP^^g^O 


but from equation (37) or (39) it can be shown that the first 
two terms in the integrand yield equal integrals and the 
minimizing condition becomes 



[P'(x,R,e)+n]Sgdx = 0 
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Since this latter equation must be satisfied by all possible 
variations of the displacement function g(x,6), it follows 
that the desired condition is 

P'(x,R,d)+n = 0 (47) 

Stated in words, the condition for minimum wave drag of a 
quasi-cylindrical body of given volume is that the longi- 
tudinal gradient of pressure on the body in the combined 
forward and reverse flow field is a constant. Furthermore, 
from equation (44), minimum drag is then given by 

D ntn =%(y-M)=%V. (48) 

where p can now be identified with the negative pressure 
gradient in the combined field and (V—irR?l)=V, is the 
volume exposed to the fluid around the cylindrical control 
surface r—R. 

The actual cross-section-area variation of a minimum- 
drag constant-volume body and its pressure distfibution 
are shown in figure 5 for the case in which axial symmetry 


r 



Figure 5. — Mimimum drag body with pressure distributions in 

combined flow. 


is imposed. Pressure coefficient in the combined flow field 
of an axially symmetric body has been given in equation 
(42). The geometric criterion just established then leads 
to the integral equation 

2S ' (i >-X' s ' ( *'K^)§-- <“+» < 49 > 


DRAG MINIMIZATION 

In this division, optimum bodies having certain prescribed 
geometric properties will be determined by standard varia- 
tional methods. The analysis will, as mentioned previously, 
deal with the strength of an axial source distribution as the 
minimizing function, rather than the geometric quantity, 


r 



Figure 6. — Body and associated nomenclature used In drag 
minimization. 


cross-sectional area. Thus, we shall be concerned with 
formula (32), giving drag in terms of the source distribution. 

QUASI-CYLINDRICAL BODY OF REVOLUTION OF GIVEN VOLUME 

Isoperimetric conditions. — The configuration to be con- 
sidered, together with associated nomenclature, is shown in 
figure 6. The geometric properties of the body can bo 
expressed in terms of the source distribution function A 0 (x) 
by using equation (20), namely, 


C, \ 2 irR 

S ( A=-Jj ~ <Pr 


_ r-R 


Then, from equation (30) 

»/.a 1 r-** aj(x i)(*-*i) 

If equation (51) is integrated x-wise it is seen that 


(50) 


(51) 


and the solution of this equation will determine the body 
geometry. In the following section an analogous integral 
equation will be derived but with the source-strength dis- 
tribution chosen as the fundamental dependent variable. 

The combined-flow-field technique can also be used to 
study the problem of minimizing wave drag for specified 
body caliber or, more generally, when the body has a fixed 
cross section at a specified longitudinal position. The 
resulting condition for minimum drag is that the pressure 
distribution on the body in the combined flow field is a 
constant forward and aft of the specified position. These 
conditions are all analogous to those obtained for planar 
problems by R. T. Jones (ref. 13). 


S(x)-S(o)=-j- 


~f> R (z—x l )A c '(x 1 )dx i 
PR V(z— Xi) S — 


(62a) 


By changing the order of integration and performing the 
integration with respect to z one finds 




x ~ pR jx—x^A^Xi) 
pr -j(x—xi y—fpR 2 


dx i 


(52b) 


or, integrating by parts, 

1 /» x—fiR 

S(x)-S(o)=±- 4/(*i)V(*“*i y-pR'dxt (52c) 

UoJ -0R 
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The magnitude of the additional volume wrapped around 
the cylinder is 

V e =j' o [S(x)-S(o)]dx (53a) 

and, from equation (52b), 

i ri-prt , 

V t =±- A 0 (x 1 m-x i y~P i R i dx 1 (53b) 

U 0 J -pR 

The variational problem. — The quantity to be minimized 
will bo taken as D—piV,. From equation (32), the drag 
can bo written, after an integration by parts, 

D= Po r-fi R MxM* C 1 -** A/faWl-XiY-pR* 

4ir J _/j b ya—x^f—^R" 1 J -j)k Xi—Xi 

(54) 


is a condition to be met by the minimizing function A 0 (x 1 ) 1 
and by its variations. 

The quantity D — mi V e can be formed from equations 
(53) and (54), and if the variation is performed, one finds the 
condition 

C'-f* SAJxQdx i f (‘'-exA o '(x i )^(J,-x 1 ) i -0 l E l , 

J-pb \J_, S zm “** 

If this last equation is compared with equation (55), it is 
seen that for admissible variations, the quantity within the 
brackets must be set equal to \(l—x i), where X is an arbitrary 
constant. Thus, the equation for determination of the 
optimizing source distribution under the conditions of given 
volume and closure is 


In addition to prescribing the volume added to the funda- 
mental cylinder, we shall also require that the body return at 
the end to the same cross-sectional area as at the front. 
Thus, according to equation (52b), 


(l—xJApjxJ dxi 

-PR V(* X lf 0‘R? 


(55) 


X 


l-fiR 

-PR 


— 7 - g - J ^ dx l =Ul-x)+^ [(l-xf-m 

(56) 


Equation (56) is recognized as the familiar airfoil equation 
with [A 0 '(x 1 )y(L—x 1 f—p 2 R 2 \ as the unknown. Thus we 
write the solution immediately as (see, e.g., ref. 15) 




j: 




~ u 1 pjj, 


X—Xi 


-^j(l~pR—x{) {xx+pRjdxx 


The first integral on the right vanishes according to the 
closure condition (55) and, if the remaining integrations are 
performed, we find 


A/wa-xy-FR 1 


i 

tttJ(1—PR—x)(x-\-PR) 


2ir] 
pJJ o 


Mi+2X 


•) 


1—fiR—x 

8 


[' 


4X(Z— 2x)+ 


2^1 
pjj 0 


(3P-4PRI -8j5*B*-12ix+8®*) 


]> 


Properties of the optimal source distribution. — It is con 
venient to express the various quantities such as source 
strength, area distribution, etc., as dimensionless functions 
of the dimensionless variable t]=x/1 and parameter <r=fiRfl. 
Thus, indicating a dimensionless function by a star, we have 
from equation (58a) 

[a-2,)V(,-l-«)(i 

^ oos " TT§] (58b > 


It will be noted that unless 

2 irl 

pjj 0 


Mi+2X=0 


this solution for A„'(x) does not obey the closure require- 
ment. Therefore we impose this last condition and finally 
obtain 

mi (P-^RlS^R^Slx+SA 


A/(x)=- 


(57) 


4pot7 0 -yJ(l-\-pR— x)(x+pR) 

The strength of the minimizing source distribution A„{x) 
is now obtained by integrating equation (57); 




(l—2x)-\J(l+pR—x)(x+pR)— 


20 , R 1 cos' 


l-2x ~| 
l+2pRJ 


(58a) 


It will be noted that if the radius of the control surface is 
taken very small, so that cr— >0, formula (58b) becomes 

Ao*(v)U «-=£r (1-2,) 

ego 

which is the well-known slender-body theory result for the 
source distribution corresponding to an optimum body of 
given volume (refs. 16 and 17.) 

In order to determine the value of the Lagrange multiplier 
Mi, in terms of the prescribed volume V e , it is convenient 
to find first the expression for the local cross-section area of 

the optimum body. Thus, using equation (52c) ^with 

S*(n)=-p Sdv))’ 
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S*(v)-S*(0)=^(r, +2a)(l- V +2a) 

[T](l — ij)E—cr(l—4:cr)(K—E)] (59) 

where K and E are elliptic integrals of the first and second 
kinds, respectively, of modulus 

>-?_ v(l—v) 

0j+2<r)(l — 17+2 <t) 

Using equation (59) in equation (53a), we find 

F «* = £r f V( 7 ?+2<r)(l— j?+2cr)[7j(l— r])E— a( 1 —4 o) (K—E)]di J 
Offo Jo 

(60a) 

which expresses the constant jun in terms of the prescribed 
volume V ,* and of a function B of the quantity 

<r=f}R/l. A graph of this function B(o) versus a is shown in 



where V 0 is the volume of the original cylinder section, 
Vo = v RH=lS(o)=PV 0 * 

Consider the expression (61) for the source function A 0 *(t]). 
In the parameter cr=f3R/l, we may think of /9 as fixed and l as 
unity, so that variations in a amount to variations in the 
size of the control surface of radius R. Thus, in figure 8, the 


4h)/v/ 



case c=0 corresponds to the source distribution for the well- 
known Sears-Haack body (refs. 16 and 17). It will be noted 
in the cases where <r> 0 that the source functions become 
less steep and attain lesser maximum values because the 
volume remains the same while the control-surface cylinder 
is increasing, thus giving a smaller maximum radius of the 
added portion. 

Next let us examine the expression (62a) for cross-sectional 
area. First, we notice that it can be written 

S*(v)-S*(o-)=^ y(i?~b2<rl(l 7]-j-2tr) 

[tj(1 — rj)E—<y(l — 4o-)(£l — E)\ (62b) 

in which form it reduces formally for a — >0 to 

<s*(i/)=^ v*{vii-v} w mo) 


figure 7. Shown also in figure 7 is a dashed line that cor- 
responds to the asymptotic form for B(cr), which is 

«')~s (60b) 


The closeness of the asymptotic values to the exact values 
even for relatively small values of a is noteworthy. 

The formulae (58b) and (59) for the source strength and 
cross-section area, respectively, can now be recast in terms of 
prescribed quantities 

[a -2v) VOrMO-Tj+v)- 2 ^ cos- 1 !=g] 

( 61 ) 


S*W-S*(o) (V*IV*) 
S*(o) B(a) 


V(v+ 2 tr)(l — >7+2 a) 


which is identical with the expression for cross-section area 
of a slender optimum body of prescribed volume (Sears- 
Haack body). Of course, V, is, in this case, the total volume 
of the body. On the other hand, if we allow the radius of 
the control surface to increase indefinitely, equation (62b) 
gives (using the asymptotic form for B(a), oq. (60b)) 

S*(v)-S*(o)=6V t * v (l-v) 

In the case when R is very large, we take 

S(x)-S(o)=2rRAr(x) (63) 

so we have, returning to the original variables, 

*<*>-<> 2 as ^ ' 64 > 


[17(1 — 17 )E — cr(l — 4 <t)(-ET — Ej] (62a) 


where V,/2 *Rl is a finite quantity, and, in fact, is the average 
height of the protuberance above the control cylinder. This 
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result is clear from physical reasoning, for one would expect 
that as the control cylinder increased in radius, the two- 
dimensional result for the optimum problem would become 
more nearly valid, and, indeed, equation (64) is the formula 
for a two-dimensional biconvex section, where Ar is dis- 
tance from the mean line, l is chord length, and maximum 

thickness is 

It will be noted that the area distribution as given by 
equation (62b) has fore-and-aft symmetry, since the function- 
al dependence upon rj involves only the combination 77(1 — u)- 
The maximum cross-section of the optimum body then oc- 
curs at the midpoint y=% and is given by (from eq. (62b)) 

S* mat—S ¥ {o)=2V* ijg [f - <7(1 ~4cr){K~E)j=2V* T(a) 

(65) 

where the modulus of the elliptic integrals is now k= 1/(1 +4 a ) . 
Figure 9 shows the function T{cr) versus a. 




Figure 9.— The function TO). 

The drag of the optimum bodies can now be evaluated. 
From equations (54) and (56) 

®-f; £7 w=§qm } 


The integral involving X vanishes because of the closure 
condition. The remaining integration gives 

D= W 0 XT Ao( ° ( 66 ) 


by equation (53b). Finally, using the evaluation of px of 
equation (60a), we have 


D o V. 2 
So - a m<T) 


(67) 


Numerical results pertaining to the problem just solved will 
be given in a later section, and a summary of the important 
formulae is given in the appendix. 


QUASI-CYLINDRICAL BODY OF REVOLUTION WITH GIVEN CALIBER 


The variational problem. — For this problem, we prescribe 
the area at the base of the body, so the given condition is, 
from equation (52b) 


AS=S(l)-S(o)~J l _ 


-f> B (l—x 1 )A 0 (x 1 )dx 1 
SR -y/(Z— r,)*— /3 2 iZ 2 


( 68 ) 


The variation can be taken as before (now without invoking 
the closure condition) on the quantity D+XAS, and it 
leads to the integral equation 

1» „-,) (69) 

J —PR % *1 PcfJ o 

The solution to equation (69) can be written immediately 
by analogy with that for equation (56). The presence of a 
linear singularity in A 0 ' (x) must be disallowed, however, 
according to a condition mentioned in deriving the drag 
formula, equation (32). Thus, we have here to set 


8g 0 A S 
orZ(Z+4jS.R) 


and the solution consistent with the given conditions is 


4 Uo(AS) l—2x 

UfHpRj Tj(l+f}R-x)(x+pK) 


(70a) 


Integrating this expression, we find for the strength of the 
optimizing source distribution 


A a (x)=l Vf^ M+fiR-xXx+fiR) (70b) 


The source distribution of equation (70b) represents the first 
approximation to the result of reference 3 for nearly equal 
front and rear radii. 

Properties of the solution. — As in the section on the body 
with prescribed volume, we now consider x made dimension- 
less by division by Z, and again set v—pR/l. The various 
quantities of interest in connection with the caliber problem 
then become 


A*(v)=l (70c) 


S*(y)-S*(0) 2 (!+2<r )(l+4<r)n(a 3 ,Z:)-(l-27;)(77+2v)(l-7 7 +2orW-(l+2cr)(l-77+2v)£: 

AS* tK 1+4<7) V(’)+2v)(l-7?+2<r) 


where n (a 2 , k) is a complete elliptic integral of third kind of 
v(l — v) t i_ —V 


modulus P= 


(ij+2<r) (1 — 77+2cr) 


and parameter <£= 


1 — 17— j— 2cr 


Again K and E are complete elliptic integrals of the first 
and second kinds, respectively, of the same modulus k. 

If we allow a to approach zero, equation (71) becomes, in 
the limit, 


|^=?[sin- 1 V^-(l-27 ? )V^(I =: ^] (72) 

which is the shape function for the well-known K&rmtin 
ogive (ref. 18). At the other limit, when a — , equation (71) 
gives (in the original variables) 

S(x)—S(o)_x 
AS ~l 
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or, using the approximation of equation (63), 

Ar(x) x 

A r(l) l 


(73) 


which is again the expected two-dimensional result for speci- 
fied caliber. 

The drag can be found by substituting equation (69) in 
equation (64), and then using equation (52b). There results 


D 4 (A S)> 

q 0 ?r(l+4<r) l 2 


(74) 


A summary of formulae pertaining to this body will be found 
in the Appendix. 


EXAMPLES OF OPTIMUM BODIES 

I 

The optimum body of given volume. — In order to examine 
in detail the dependence of the body geometry on the param- 
eter cr, we may return to equation (62a). The quantity 
S*(o)] is actually the local cross-sectional area added 
to the basic cylinder by the action of the source distribution. 
In figure 10 are shown some cases of optimum bodies, having 
equal additional volume V e *, for several values of the param- 
eter cr. Only half of each distribution is shown, since they 
are symmetric about the point 17 = 1 / 2 . The one labeled 
o=0 is the Sears-Haack optimum body, and it will be noted 
that as a increases, the curves depart rather quickly from 
this limiting case and approach the other limiting value of 
the biconvex distribution for «r-> «. In fact, a biconvex 
arc drawn through the end points of the o=l /2 case is indis- 
tinguishable from the exact result in the scale used. In the 
inset of figure 10 is shown the variation of the drag of the 
optimum bodies as a function of a. This drag is also based 
on equal volume, and shows a fairly rapid decrease with 
increasing values of cr, due to the decrease in the thickness 
of the exposed portion of the body. The dashed curve on 

the drag plot is the calculated drag 

the assumption that each meridian section of the body acts 
as an independent two-dimensional optimum airfoil. This 
admittedly crude approximation is of course very poor at 
low values of cr, but its accuracy becomes surprisingly good 
for cr greater than about 0.4, and the approximation becomes 
exact in the limit <r—>°° . 

The variation of local cross section with 0 - can be examined 
also on the basis of equal exposed area. Thus, using equa- 
tion (65) in combination with equation (62b), we have 

■ AS\bu = 2T(<t)B ( 0 ) V(’ 7 + 2 <r)(l-T 7 + 2 cr) 

[17(1 — t])E—(t(1 — 4cr) (K—E)] (75) 


D 12 , 

————— — under 

UoWZ 2 ) 2 rrcrj 


Figure 11 shows plots of equation (75), and it is again noted 
that the departure from the slender-body approximation 
(cr=0) is rapid. The limiting variation of area for cr->«> is 
also shown in figure 11 , and it is seen again how closely the 
optimum body-shape functions approach this limiting result 
even for moderate values of cr. Also shown is the drag 
corresponding to these cases. 


D 3 

4B(a)[T(a)Y 


(76) 


which shows a similar drop from the 0 = 0 value as cr is 
increased. Again, the effective fineness ratio of the bodies 
is increasing with cr, and, if frontal area exposed to the 
stream is held fixed, the maximum thickness of the excres- 
cence vanishes as 1/c r for large cr. The departure of the 
geometric variation from the slender-body case is most pro- 
nounced near the nose, i?= 0 , where the slope is given by 


dS*(y) 3tt V* 
dt] 8 B(a) 



(77) 


which vanishes only qs Vo- for er-» 0. 

The optimum body of given caliber. — In this case, the maxi- 
mum cross-section occurs at 77=1 so there is no longitudinal 
symmetry. Figure 12 shows, for several values of the pa- 
rameter cr, the optimum, equal-caliber, incremental cross- 
section area given by equation (71). The inset shows the 
drag as a function of cr; from equation (74) 

Again in this case, the closeness of the optimum distributions 
as cr increases to the two-dimensional value (a — >°°) is 
noticeable. This point has also been made by Ferrari in 
reference 5 where problems similar to ours are treated 
by a different approach. If the expression for cross- 
section area (eq. (71)) is expanded in powers of 1/cr, it is 
found that 


AS*(v,T) 77 ( 1 -^( 1 - 277 ) AY, 
AS*(l,<r) V 32 W 


(78) 


which shows the smallness of the correction to the two- 
dimensional result for moderate values of cr. 


BECIPBOCITY BELATIONS 

The optimum body of given volume. — The longitudinal and 
radial perturbation velocities can be determined by substi- 
tuting the derivative of the source-distribution function (oq. 
(57)) into the formulae (29) and (30). We find, at any 
point (x, r ) ( r>R ) in the field 


Vx (g£) _ Mi 1 

Uo 8irq 0 Vb+ZKr+A)] [Z-x+j3(r+I?)] 


{4[x+p(r+R)] [l- x+ p(r+R)\E+ 


[l(l+4£K)-±(l+2pR)(l-x+pr+pR)]K+4:(l+2pR)(l-2x)TL(c?,k)} 


(79a) 
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Fiqtob 11. — Geometry and drag characteristics of optimum bodies of given volume (bodies having equal additional frontal area). 
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<Pr(x,r) _ III 1 

Uo 87T&T ^[ x+ ^ r+ E)][i~x+p(r+R)} 

{(l-2x) [x+ftr+R)] [ l-x+p{r+B)]E - 
flrl[l—2x-}-2p(r—R)]K--4:P 2 R(r~R)[l—x J cfi(r-{-R)]K-{- 
4f} i (P—R^(l+2f}K)TL(a! > ,k ) } (79b) 

where now 

x (l — 

¥== [x+p(r+R))[l-x+p(r+R)] 

X—P(r—R) 

l-x-\rP{r-\-R) 

For the present axis system, the act of reversing the flow 
amounts to substituting l—x for x, and, for the case of the 
symmetric body, the longitudinal perturbation velocity in 
the reversed flow is 


v{x,r)=—u(l—x,r) — —<p x {l—x,r) 

Now, from equation (40), pressure in the combined field is 
given by 

P= — pJJ o{u+u) = — pJJ 0 [<pJ.x,r) — <p x (l—x,r)] (80) 

In order to verify this relation using equation (79a), the 
following result, which can be derived from formulae of 
reference 20, section 400, will be useful: Let 


then 


This, together with the fact that the modulus k is invariant 
to the substitution of l—x for x, enables us to 'write im- 
mediately 

(l-2x) (81) 


Differentiating equation (81) we find 


P'+Hi=0 


(82) 


which agrees with the criterion for mi ni mum drag with 
given volume established in equation (47). The Lagrange 
multiplier m is therefore identified as the pressure gradient 
in the combined flow field. It will be noted that equations 
(81) and (82) hold everywhere within the enveloping forward 
and rearward Mach cones of the quasi-cylindrical body 
(see fig. 4). 

Now considering the radial component of perturbation 
velocity <&■, we find, using the relations mentioned just 
previously 


Vrfar) , <pr(l—x,r) 0 \r 3 —R 2 ) 
U 0 + U 0 ~ 4 qjr w 

so that the relation 


( 83 ) 


T r 1 (r)=n(o 2 ,t) 


<p T {x,r)=—<p T (l—x,r) 
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is satisfied on the quasi-cylinder itself, that is, when we set 

r=R. 

The optimum body of given caliber. — For this case we find 
the following equations for the perturbation velocities: 

Vx( x i r ) \K— OTTfVv* 1A 1 

U. ~ 2^[x+f)(r+R.)}ll-x+fKr+R)] L 1 M 

(84a) 

<Pr(x,r) X 

Uo 2icq<?^[x-\-fi(r-\~R)] [l—x+&{r+R)\ 

{ [x J t-p(r+R)][l—x+p(r+R)\E—pr{l+2pR)K} (84b) 


where the elliptic integrals have the same modulus and 
parameter as in the previous section. 

In this case, the pressure in the combined flow field is 

P=-pJJ 0 [<p I (x,r)-\-<p x (l—x,rj\ (85) 

which gives 

P = — X (86) 

so that in this instance, pressure itself is constant in the 

combined flow field. 

From equation (84b) , we see that 

Vr (a:,r)= ¥ v(Z— x,r) (87) 


since the modulus of the elliptic integrals is invariant to the 
change x-^l—x. 

Uses of reciprocity relations. — The reciprocity relations 
serve the dual function of checking the derived perturbation 
potential against minimiz ation criteria based on other 
considerations (see eq. (47)) and of relating the Lagrangian 
multipliers to the pressure or pressure gradient in the 
combined flow field. Equations (81) and (86) also reveal 
that the expressions for pressure in the combined flow field 
hold, independently of r, throughout the entire region within 
the enveloping cones of the bodies. These results are 
generalizations of a similar effect noted in reference 19, 
where the combined pressure field associated with a Sears- 
Haack body was shown to have a constant gradient within 
the enveloping cones. In the latter reference, this property 
of the minimum-drag bddy was used to expedite the calcula- 
tion of interference drag with a satellite body lying within 
the enveloping cones. Similar methods could obviously be 
applied to the present configurations. 
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APPENDIX A 


SUMMARY OP FORMULAE FOR THE OPTIMUM BODIES 


The formulae derived in the text for the body shape func- 
tion, pressure coefficient, and drag of optimum bodies having 
given volume or given caliber are repeated here for con- 
venience. The type of configuration treated, and the nomen- 
clature, are shown in figure 13. 


r 



Figure 13. — Body and associated nomenclature used in drag 
minimization. 


THE OPTIMUM BODY OF GIVEN VOLUME 

The variation of AS for the optimum body with given 
volume is 

NS(x) = V / R p\ V (g+2 PR) (l—z+2pR) 

“*(?) 

{x(l-z)E(/c)-pR(l-4m [K(k) —E(k )] } (Al) 

where 

A S(x) =Tr[(R+Ar) a -R?\ 

V, = volume of exposed portion 

function defined in equation (60a) and shown 
in figure 7 

K(k) =complete elliptic integral of first kind of modulus k 
E(k) = complete elliptic integral of second kind of modulus 
k 

u _ x(l—x) 

(x+2pR)(l— X+20R) 

Examples of optimum bodies for a few values of the param- 
eter f}R/l are shown in figures 10 and 11. 


The pressure coefficient on the body is 



u _ V t v 4(x+2pR)(l-x+2pR)E-[l(l+4pR)-Hl+2pR)(l-x+2pR)]K+4:(l+2l3R)(.l-2xmcf,k ) 

U° 2tt lKB (Mj ■yl(x+2pR)(l-x+2pR) 
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where IT (a 2 ,/:) is a complete elliptic integral of third kind of 
modulus k and parameter a 2 (in the notation of ref. 20). The 
parameter a 3 is given by 


1-X+20R 

Figure 14 shows some plots of C P /(VJP) versus x/l for a few 
values of the parameter pR/l. 

The wave drag of this optimum body is given by 



(A3) 


The variation of drag with /3 R/l is shown in figures 10 and 11. 

THE OPTIMUM BODY OF GIVEN CALIBER 

The variation of A S(x) for the op timum body of given 
caliber is 



body of given volume. 


A a/_-. 2 A S(l) ^J(l+2pR)(l+40R)IL(ci*,k)-(l-2x)(x+2pR)(l-x+2pR)E-l(l+2pR)(l-x+2pR)K 
{) *i(i+ mr ^(x+m)(i- x +2W) (A4 > 


where the symbols have been defined above. Examples of 
optimum bodies for a few values of the parameter /S R/l are 
shown in figure 12. 

The pressure coefficient on the body is given by 


, = _S A S(l) 1+20R 2n(c?,k)-K 

P l l+mt -yJ(x++^R)(r—x+2(iR) 


(A5) 



Figure 15. — Pressure distributions for various cases of the optimum 
body of given caliber. 


Figure 15 shows some plots of 

values of the parameter /3 R/l. 
The drag of this body is 


versus x/l for several 


D 4 [A S(k)] 3 

?0 v(l+ / if}R) l 

and its variation with J3 R/l is shown in figure 12. • 


(A6) 
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